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Abstract 

Inspired by the recent HEIDELBERG-MOSCOW double beta decay ex- 
periment, we discuss the neutrinoless double beta decay in the super- 
symmetric seesaw model. Our numerical analysis indicates that we can 
naturally explain the data of the observed neutrinoless double beta de- 
cay, as well as that of the solar and atmospheric neutrino experiments 
with at least one Majorana-like sneutrino of middle energy scale in the 
model. 
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1 Introduction 



It is commonly considered[l] that neutrinoless double beta decay ((/3/3)oi/) is a very sensitive probe for 
new physics beyond the standard model (SM). The models, such as the Majorana mass of light neutrinos[2], 
right-handed weak couplings involving heavy Majorana neutrinos[3, 4], the Higgs-boson contribution to {(3(3)qu 
in an SU(2) x U(l) gauge model with left-handed Majorana neutrinos[5], as well as the R-parity violation 
supersymmetric model that was first proposed by Mohapatra[6] and later studied in detail by other groups[7, 8, 
9, 10], can result in the observed neutrinoless double beta decay. Definitely, the steadily improved experimental 
bounds[ll] on the {(3(3)qu lifetime can then be translated into more stringent limits[12] for the parameter space 
of these models. This would be extremely valuable information to our search for new physics beyond SM. 

Recently, a positive indication of the neutrinoless double beta decay has been reported[13]. According to 
the announcement of the Heidelberg group, the half-hfe time of {l3(3)oy for nuclei ^*^Ge is (0.8 - 18.3) x lO^^ 
years (with the best value of 1.5 x 10^^). Many Physicists focus their attention on the work due to its 
important consequences to the particle physics and astrophysics. Barger et al. [14] pointed out that accurate 
measurements of neutrinoless double beta decay may constrain the neutrino component of the dark matter. 
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Combining the bound set by the CHOOZ reactor, authors of Ref.[15] derived constraints on the neutrino 
mixing angles; a model-independent constraint on the neutrino mass spectrum imposed by the HEIDELBERG- 
MOSCOW double beta decay experiment is also given in Ref.[16]. Assuming a reasonable hierarchy pattern of 
neutrino masses, the magnitude of neutrinoless double beta decay is estimated in Ref.[17]; Prigerio, Smirnov 
and other authors [18] evaluated the effects of various neutrino-mass-matrix structures on the neutrinoless 
double beta decay. The authors of Ref.[19, 20, 21, 22, 23, 24] discussed neutrinoless double beta decay in 
different possible models. In Ref. [25] , the effects of the leptonic CP violation phase on the neutrinoless double 
decay have been investigated. Meanwhile, discussions on the neutrinoless double decay in nuclei ^^^X^ are 
presented in Ref. [26]. 

By contraries, several groups suspect if the present data can definitely indicate a non-zero rate of (PP)ou- 
Aalseth et al. [27] pointed out that extraction of those signals depends on the choice of window, and some 
preset conditions, such as the absence of a flat background and the relative strength of the ^^^Bi peaks etc. 
The authors of Ref. [28] also criticize the claim of "evidence" which depends on the data-set choice. The most 
important point is that a previous analysis[29] does not find any hint towards (PP)ou with the same data. 
Even so, we are inclined to believe that the conclusion about non-zero rate of (/3/3)oi/ is positive and then we 
can extract constraint on the parameter space of "new physics" models from the obtained experimental data. 

As it turns out [30], we cannot explain the solar neutrino, atmospheric neutrino together with LSND 
collaboration observes with three flavors of neutrinos. In the literature, the puzzle is addressed through the 
following approaches. First, one of the three observations is simply discarded. In the second approach a 
"sterile" neutrino[31] is introduced and its existence does not affect the decay width of Z boson. In the 
last approach there exists a mass-discrepancy between neutrino and anti-neutrino which is realized by the 
CPT violation terms in the neutrino sector [32]. Following a more common point (which is set without any 
justification), we adopt the first approach while ignoring the LSND observation due to its relatively large 
experimental error (more than 3 a) here. 

If we believe the announcement of the non-zero rate of (/3/3)oi/ and assume the light-neutrino-exchange is 
the dominant mechanism for the process, by considering matrix elements of Ref. [33], which include the con- 
tributions from higher order terms of the nucleon current, one finds that to fit the HEIDELBERG-MOSCOW 
experimental results the effective neutrino mass {rriu) = (0.11 — 0.56)eV (95% c.l.), with the best value of 
0.39cV. Though there are a few neutrino mixing schemes which can simultaneously explain the data of the 
solar and atmospheric neutrino experiments within the three generation model[34], the most favorable sce- 
narios of neutrino mixing lead to {nii,) < 0.01eV[35]. This discrepancy implies that another mechanism might 
be responsible for the neutrinoless double beta decay process. 

In this work, wc will discuss the loop-induced (/3/3)oi/ in the minimal supersymmetric extension of the 
standard model with right-handed neutrinos (MSSMRN). Although the seesaw mechanism may lead to the 
non-zero Majorana masses for three light neutrinos [36], the contributions from the virtual Majorana neutrino- 
mediate diagrams are suppressed by a small factor -^^^[37], where (mj,) is the "effective Majorana mass 
parameter": |(m,y)| = lU^imi^^ +11^21^1^^ +U^^mi,^\ and m^. denotes the mass of the Majorana neutrino Uj and 
Uej the clement of the neutrino mixing matrix, pp ~ lOOMeV is the nucleon Fermi momentum. If assuming 
the constraint on < > as (m^) < O.OleV (90% C.L.) [38], we find ^ < 10"^°. On the other hand, the 

Pf 

loop diagram contributions may play an important role in the neutrinoless double decay if we have a relatively 
light Majorana-like sneutrino. A point should be noted that a similar computation has been performed by 
Hirsch et al.[39] in the supersymmetric model with non-universal soft breaking terms and mass terms of the 
Majorana- type neutrino. The model we adopted here is a concrete realization of the general case discussed 
in Ref. [39]. 

The paper is organized as follows. In section 2, we review the minimal supersymmetric extension of the 
standard model with right-handed neutrinos and give the notations adopted in our analysis. In Section 3 we 
derive the supersymmetric contributions to the AL = 2 effective lagrangian at the quark level. The derivation 
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of [(3(3)^^ transition operators and nuclear matrix elements are given in section 4. Under some assumptions 

for the supersymmetric parameters, we give the numerical analysis on ( ) decay of ^^Gg in section 5. Our 
conclusions and discussions are made in section 6. Some complicate and tedious formulas are collected in the 
appendices. 



2 The supersymmetric extension of the standard model with right- 
handed neutrinos (MSSMRN) 

As the right-handed neutrinos are introduced into the game, the superpotential is written as[40] 

+h1^e,iH}Q]D' . (1) 

Here, ff^, are the Higgs superfields, L^, are the superfields in doublets of the weak SU(2) group, where 
1=1, 2, 3 are the indices of generations, the rest superfields , , and W are in singlets of the weak 
SU(2). Indices i, j are contracted for the SU(2) group, and /i", h'^, /l^ h" are the Yukawa coupling constants. 
To break supersymmetry, non-universal soft breaking terms are introduced as 



+AleijHfQ'^U'^ + h.c], (2) 
where , , m? „, rn? rn^r,-, fn^r,, rn^r,^ "7,^,, and m^,, are parameters in unit of mass square, 

/fl ff2 ' ifl2 ijIJ ijIJ j^IJ jiU i^IJ ^ -L J 

ma, 7742, JT^i denote the masses of (a = 1, 2, 8), \\ [i = 1, 2, 3) and Ab, which 

arc the S'C/(3) x SU{2) x U{1) gauginos respectively. B" are free parameters in unit of mass square. 
Ajj, Ajj Ajj, Ajj (/, J = 1, 2, 3) are the soft breaking parameters that result in mass splitting be- 
tween standard particles and their supersymmetric partners. Taking into account the soft breaking terms 
in Eq.(2), we can study the phenomenology within the minimal supersymmetric extension of the standard 
model with right-handed neutrinos(MSSMRN). The gauge symmetry SU{2) x ?7y(l) breaks down into U^{1) 
through the nonzero vacua of two Higgs fields. Different from the MSSM, neutrinos of three generations 
obtain nonzero Majorana masses through the seesaw mechanism [40] 

(hi \ ( K \ 

/i^ (m-)-i h'i U, (3) 

V hi) \ hi) 

where the right-handed neutrino mass matrix is introduced in Eq.l and the unitary matrix U is used to 
diagonalize the neutrino Yukawa coupling matrix h^^ . 

In general, the mass matrix m,^ is still not diagonal in the weak basis, so that we need another unitary 
matrix to diagonalize ihjy 
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(4) 



Thus as far as UU^ is not trivially a unity matrix I, the mass eigenvalues are non-degenerate, one can expect 
neutrino oscillation among different flavors which is a target of current and future experiments. The mixing 
for L{ and , N^* results in nine scalar neutrinos. With the basis = A''^, N^*), the resultant mass 
matrix of the scalar neutrinos is written as 



with 



Pi 



Pij 

\{rrf'^m'')ij + hfvl5, 



IJ 



Pi 



IJ 



\ 



(5) 
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C^2 



C2 
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(6) 



where = cosO^^, = sin^^ and 9^^ is the Weinberg angle; vi, V2 are the nonzero vacuum expectation 
values of the two Higgs sectors and tan /3 = ^ is the ratio of the two vacuum expectation values (hereafter 
we employ the symbols = cos/3, = sin/3). The eigenvectors L{, , N^* of the weak interaction form 
nine sneutrino mass eigenvectors P {i = 1, 2, • • • , 9) by diagonalizing the matrix m|: 



= {ZlfL{ + (Zt)^(3+-f)iv-f + (4)^(6+^) iV^* , 



(7) 



where the L{ are the Dirac-type sfermions and the Lagrangian containing L{ conserves the lepton number, 
whereas that containing the Majorana-like scalar violates the lepton number. As for the charged sleptons, 
fields L2 and mix to produce six mass eigenvectors {i = 1, 2, • • • , 6) of the squared- mass matrix: 



/ 



\ 



^ {hjiJ,V2Sij + A.fjVi^ 



\ 



mz 



Zl = diag(m^__ , 



For the gaugino and scalar quark sectors, we take the notations of Ref.[41]. Now let us turn to the 
phenomenology of the super symmetric model with right-handed neutrinos. 



3 {/3^)qj^ decay in MSSMRN: the effective Lagrangian at quark level 

In this section, considering the contributions of the model discussed above to ^ , we derive the relevant 
effective Lagrangian in terms of color-singlet currents. 
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Within the framework of MSSMRN, there are two sources of lepton number violation, the contributions 
from the Majorana neutrinos mediated processes and the Majorana-hke sneutrinos mediated processes to the 
( PP) [8]. At the quark level, (fSfS) decay is induced by the transition of two d— quarks into two u— quarks 

plus two electrons. All possible Feynman diagrams that induce the (fSfS) decay are shown in Fig.l. As 

stated in the introduction, the contributions of Fig. 1(a) to (^P/^j^ decay is suppressed by a factor 

Besides this suppression factor, the contributions of Fig.l(b,c,d) to {(3/3) still suffer from the typical loop 

V / Ou 

integration suppression. In contrast to the SM case, the supersymmetric contributions are different. Three 
classes of Feynman diagrams are given in Fig.l(e,f,g). It is noted that the supersymmetric contributions only 
receive the loop integration suppression which indeed is unavoidable. Therefore, we can ignore the SM sector 
in our later calculations, but only keep the supersymmetric contributions, namely in the following section, 

we only consider the three classes of contributions (Fig.l(e,f,g)) to the {(3f3) decay. The explicit Feynman 

V / Ou 

diagrams of the three classes at one-loop level are drawn in Fig.2 and Fig.3. In those figures, Fig.2(a,b) belong 
to the class shown in Fig.l(e); Fig.2(c,d,e) and Fig.3(h,i) belong to the class shown in Fig.l(f) and Fig.2(f,g) 
and Fig.3(j,k,l,m) belong to the class shown in Fig. 1(g). 

After integrating out those heavy internal particles, the effective Lagrangian for dd — uu+e~e~ is written 

as 



''ff ~ 2mp\ 



+C2(^w)^a7p(l - 75)C?Q%7''(1 - 75)1^/36(1 - 75)6*= 

+C3{iJ.^)ua{l + 75)c?a^i/37p(l " l5)dpeY{l - 75)e'' 
+C4(/x^)na(l - 75)c^a^t/37p(l " l5)df)eY{l - 75)e'' 
+C5(/Ww)^a(l + 75)c^a%7p(l " l5)dpeY{^ + 75)e'' 

+C7{iJ,^)(uaC0+daUpu!+dp + ^ua ^yUJ+daUj3a^''' LO+dp^euj-e'' 

+ Cs{Hy^)ualp{l - J5)daupja{^ - l5)dpe{iaf^){l + 75)6^^ 
+C9(/^w)^a7p(l - 75)daUi3'Ja{^ - 75)d/3e(icr'"^) (1 - 75)6'' 

+Cio{lJ,yv)ua{i(Tpa)uj+daU/3Y^-d/3ej''u+e'' 
+Cu{fi^)ua{iapa)LU-.daUpY^-df3e'y''uj-e'' 
+Ci2{fJ-^)ua{icrpa)LO+daU/3Y^-d/3e'y'^iO-e'' 

{fJ'w)^iJ.pTSUal'^i^-daUf3crP''u;+d(3e'y^u;+e'' 

+Ci3{n^)ef,paSUa7''oj-daU/3a^uj-d/3ej^u;-e'' 
+Cii{li^)e^p„SUal^ijj-daUfjaf"LO+dfie^^uj-e'' 

-lveiilJ-w)^iJ.P'TSy'al'^^^-daUpY^-di3ea'''^uj-e'' 
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(9) 



The explicit expression of rii{fi^) and Cj{^^) are given in appendix B. We integrate out those heavy fields 
at the electro- weak scale ji^. Here, we have used the Fierz transformation to re-arrange the quark currents 
into color-singlet forms. 

The next step of the calculations is to re-formulate the transition matrix amplitudes where all quark 
degrees of freedom are replaced by the corresponding nucleon degrees of freedom in nuclei and it is based on 
the principle so-called as " quark-hadron duality". To carry out these calculations, one needs concrete forms 
of the relevant nuclear structure, i.e. the wavefunctions which describe the nuclei. 



4 Nuclear {(3/3)^^^ decay 

At the nuclear level, the amplitude for (^PP^^ is written as 

{{A,Z + 2),2e-\S -1\{A,Z)) = {{A,Z + 2),2e~\Texp\i f d'^xC' 



ff 



(10) 



where the effective Cf^^"^^'^ is given in Eq.(9). The nuclear structure is involved via the initial {A,Z) and 
the final (A, Z + 2) nuclear states which have the same atomic weight A, but different electric charges Z 
and Z + 2. The standard framework for the calculation of the nuclear matrix elements is the nonrelativistic 
impulse approximation (NRIA) [2] . Now, we turn to the derivation of the nuclear matrix elements of the color 
singlet currents in Eq. (9) using the results of Ref . [42] . The relevant matrix elements of those currents are 

{Pip)\ud\N{p')) = FP{q^)M{p)nM{p') , 
{P{p)\u^5d\N{p')) = FP{q^)N{p)^^T+N{p') , 

{P{p)\uY{l - l^)d\N{p')) = M{p)(Fv{q^) - FAiq^H + Fw {q^)icjp'' q„ + Fp{q^)qp^^)T+M{p') , 
{P{p)\uYil + l^)d\N{p')) = N{p){Fv{q^) + FAiq^H + Fw{q^)iaf\^ - Fp{q^)qp^^)T+N{p') , 
(P(p)lna^'^(l + 75)|iV(p')) =J^{P){J'"' + '-eP'^^'JxsY+Nip') , 
{P{p)\ua'^{l - i^)\N{p')) =M{p)[j'^ - '-eP'^^'jxsy^ip') , 



(11) 



where M ■ 



N 



is a nucleon isodoublet, and q = p — p' ■ The tensor structure is defined as 



mp 



For all form factors, we take the dipole form[43] 



fv,A Fs,P,w{^) T^^\0) 



m. 



(12) 



(13) 
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with rriA = 0.85GeV and /y 1, /a ~ 1.261. The nuclear couphngs obey the relations 

FwiO) _ Hp- Hn ^ 3.7 

r^w 

fv 2mp 2mp 

Fp(0) _ 2m p 
fA ml 



(14) 



where mp^^^-^ is the proton (pion) mass and is the proton (neutron) magnetic moment. With the 

nonrelativistic quark model, t[^\q) = 1.45, T^^ = -1.48, = -0.66 [42]. 
The transition operators are defined as 



a+b 

\rab ' 



aj^b 



^T' = ^a^bl^i^o. ■ fa&)(0-6 • hb) - <^a ' (^b] ( — ^F5{xa) 
a^b ^^"-^^ 

^V,N = Y^t:^T^ CFbl — ^FNixA) , 
a^b 

^V,l =YT^T'^<^f, X rab( — ^F6{xA) , 



a^b 



Tab J 



^v,2 = Y.-trt<ya[^)F.{xA) 

a^b ^^"^^ 



n 



V,3 



Y:r^r^{3{aa-rab)rab}(^)F,{xA) • (15) 

a^b ^^"^^ 



Here, Rq is the nuclear radius being introduced to make the matrix elements dimensionless and other notations 
are: 

ra6 = (ra-r6), rab = \rab\, ^ab = —, XA = mArab, (16) 

Tab 

where is the coordinate of the f-th nucleon. The above matrix elements have been written in the closure 
approximation, which is well satisfied due to the large masses of the inter-loop particles. The expressions for 
those structure functions Fi are given [39] 

i^Jv(a;) = ^(3 + 3x + x2)e-^ 
F4{x) = ^{3 + 3x-x^)e-- , 
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^6W = 7^(l + a^)e~^ 



With these form factors, the reaction matrix element is obtained as 



where 



7^' 



0+^0+ 



Ci{fi^)nv-A e(l+75)e^ +C2(Mw)^^v-a e(l - 75)6' 



+C3{fi^)n^, [e7o(l - 75)e'] + C^MQ^ [e7o(l - 75)e' 
+C5(Mw)^? [e7o(l + 75)e''l + C6(Mw)^5T [e(l + 75)6'= 
+C7{iJ.^)nST e(l-75)e 



/. 



^ Ff ( 
I2\mpj /a 



+ 9- 



,pk 



^V2 ^ ^Vs 



2mp /i "1 12Vmp; /i 
0^ = J^?(Ff (O)^-Ff (0)), 

= 'I^[(!r\\F\nF,M\I) - {F\^gt,n\I)\ > 

mpf/ Ff (0) x2 r^ffM^'/FlO \n^rw\^(TxlM\' 



+ 



'T!f\Q)^2 T^'\'^)T!i'\{))^ /mA\2 rl /T^ )(0)x2 2 ^ (0)rf ' (0) 



/3 



\m,iD/ Lh V f 4 / 



<mp/ Lev /a / 3 

4 rf (o)rf (0) _ l/if^\2n ^ /mAx2ri / rf (o) n2 1 rf )(o)rf (0) 

3 f\ 12 V /a J Vmp/ I- 12 V ; 



/I 



/I 



2rf^(0)T3^"^(0) 1 (F';>{fd) 



.(3) 



+ 3 



/I ^12V fA 



(3), 



The 0+ 0+ decay rate dvf^fr^'" for the process iVj(A, Z -2)^ Nf{A, Z) + a + 62 is written 



as 



>0+^0+|2 



spin 
-14 



piP2eie25(ei + £2 + -E/ - Mj)(iei(ie2(i(pi • P2) 



with 



(mo. 



(21) 



+ (C6(/i^) - C7(/x^))j^5Tf + ^|/i(ei)/i(e2)f |(Ci(/x^) - C2(/x^))0y-A 

-{c^{^i^) + C5{fi^))n^, - C4(mw)^^2 + (cel/^^) - C7(/x^))j^5t|' , 

= 2/i(ei)/i(e2)5-i(ei)5-i(e2){|(C3(/Xw) + C5(//w))^^2 + C'4(Mw)^^2f 
/i(ei)/i(e2)|'| - C2(Mw))^^v-a - {CsifiJ + C5(/x^))n? 

-C4(aIw)0^ + (C6(Mw) - <^7(Mw))^STf , 
where the expressions of /, g are taken from Ref.[2]. 

5 The input parameters and numerical analysis 

In this section, we present our numerical analysis on the neutrinoless double beta decay in the supersymmetric 
seesaw models. In order to simplify our discussion, we assume that the Yukawa interaction and soft-breaking 
part of the Lagrangian all are flavor-conserved, the flavor changing interactions are mediated by the CKM 
entries of the quark and lepton sectors. This is just the so called 'minimal flavor violation' scenario in the 
supersymmetric models. Under the assumption, there are four couplings, gi, g2, gs and the right-handed 
neutrino masses m", the Yukawa couplings for neutrinos /i^ (/ = 1, • • • , N^) and 6 -|- llAT^ free independent 
parameters for supersymmetry breaking, all together lO+lSiVg parameters to be fixed besides the lepton CKM 
matrix elements. Among them, 11 parameters for each generation appear only in the sfermion mass matrix and 
the mixing. Therefore, it is more convenient to choose eight physical masses and three mixing angles for the 
three charged sfermions in each generation as input parameters. As for the right-handed neutrino parameters, 
we also set = m2 = ■ ■ ■ = rn^^ = and the matrix = I for simplicity. The leptonic sector CKM 
matrix simply is U in this case. For each generation, we can specify the relevant input parameters as three 
scalar quark masses fn^jj, "^j^/) ^b'^ mixing angles O^j-, O^j for the squark sector, two-charged-scalar 



lepton masses and one mixing angle m.j, m.j, 9-^, three sneutrino masses 



m 



im_j for the slepton 



sector. As well, the Yukawa couplings of neutrinos h'j are also input parameters. Another scalar down-type 
quark mass is obtained through the relation 



cos 



9-rm ,-f-sin ^.,m — m = cos 6 .rin ,-|-sin 6 .,m — m + m,„ cos 23 . (22) 
u{ u[ ul u' d[ b{ di ^ ^ ' 

Assuming the relations m_j <^ rri j < m_j hold among the three scalar neutrino masses, the sneutrino mixing 



matrix is written as (accurate to order 0{{:^^)'^)) 

/ 

Zr.i = 

V 



1 


V2p, 





Pi 




1 






V2 


"2 

Pi 


1 


1 


"2 


V2 


V2 



(23) 



9 



where 



with = cos^O.jW? + s\T? 6 -j'w?' — m? — vn? iCl — Si). A point needs to be specified that Wf (I = 

1, 2, • • • , A'^q) are negative when we derived Eq.23. For the neutrinoless double /3 decay, the lepton number 
is violated, thus only the Majorana component of the neutrino or scalar neutrino is responsible for the 
process. Therefore the amplitude must be proportional to the mixing entry. Diagonalizing the mass-square 

matrix, there are three eigen-values, which are = f, ^ , m? = i(m^, + 2Wf Vn + m? , -|- -B" 



1 Ol 2 



m 

"3 



^/ = i (^'^j + '^h'fv2 + — B"^ , respectively and here m^^ + 2h}fv2 + m^^ is a positive constant. The 

corresponding eigenvector of m? is 

^1 

where L{, and iV^* are one Dirac component and two Majorana components in the mass eigenvector. The 
mixing entry is proportional to — 2^2^ — 5 ^ • If is negative, obviously the mixing is larger (in this 

case m-i < m,-i < m,~i) and the total amplitude is proportional to (— , J^l' — — )^ ~ — si. Otherwise 

„i I ^ mI ' "1 



if Bj is positive the amplitude is still proportional to but the relation becomes m-i < m-i < m^/ 

"2 13 2 

instead, and the mixing would be relatively small. To meet the data, the mixing cannot be too small. In the 
three generation neutrino case, the contribution is suppressed by the factor < rrii, > /pp ~ 10"^*^ in addition 
to the electro-weak coupling factor, so that cannot be substantial to explain the double beta decay data. In 

m - J 

this SUSY model with right-handed neutrinos, the mixing factor is if^, typically, m^^ is about several TeV, 

so if is of a medium energy scale, (about 10^ GeV), its contribution can meet the observed data on the 
double beta decay. 

Concerning the remaining relevant parameters, ^, mi, m2, m^, it is customary to use, in place of 

/i, 1712, m-i, the two chargino masses m ^ , gluino mass m- and mi. From Eq.l and Eq.2, one can 

'"1,2 " 

easily express the original parameters appearing in the Lagrangian in terms of the physical input pa- 
rameters (the physical masses and mixing angles of sfermions). For the charged sfermion sector, the 
original parameters are expressed in terms of the physical input parameters and the expressions can be 
found in Ref.[41]. We present the explicit expressions of the original parameters of the sneutrino sec- 
tor in terms of the physical input parameters in appendix A. The consequent new relevant physical pa- 
rameters are: tan 3, mi, m., m , , m.,, m.,, m.,, 9-,, 6-,, m-,, m.j, 9-,, m m,, m j, 

hq, m^, Usxs (I = 1, 2, 3) plus the SM input parameters. 

As for the SM parameters, we take = 1.78GeV, = 5GeV, mt = 174GeV, mi = 91.18GeV, mw = 
80.33GeV, ae(mw) = j^g, Cis{my^) = 0.12 at the weak scale. Choosing the mass of the right-handed neutrino 
masses is a bit tricky. In order to explain the data of the solar and atmospheric neutrino experiments, the 
effective neutrino Majorana mass < O.OleV is required in the most favored neutrino mixing scenarios. 
When setting m^, ~ 0.01 eV, if assuming the Yukawa coupling of neutrinos h^, is approximately equal to unity, 
the seesaw mechanism requires the mass of the right-handed neutrino to be m^ ~ 10^^ GeV, if ~ 0.1, we 
have ~ 10^^ GeV. For we take three right-handed neutrinos with the same mass m^, we assume that the 
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neutrino Yukawa couplings have the relation S> hi 2 in order to fit the data of the atmospheric neutrino 
and solar neutrino experiments. 

In our calculations, we take the input supersymmetric physical parameters as following 



TT 



rri- 



m. 



m., = 2TeV, 6., = 9., = -, 

d} ' (/I iji 9 ' 



m , 



m 



m . 



m. 



m - 



m 



2TeV, m. 



200GeV, 9 



m,2 = 2TeV, m.^ 



TT 

4 ""Ef ~ """^^ • ' " El ~ 2 ^ ^' 

m 



, = lO'^GeV, m , = m „ = m , = lO^^GeV, 

2 ^3 "I ^3 



/jf^ = 0, /i^ = , mi = m- = SOOGeV, m 
10 



: 200GeV, m , = 500GeV. 
For the mixing matrix of the lepton sector, we consider two possibilities, they are given as 



Ui 



/ 0.91 0.35 0.24 \ 

-0.42 0.72 0.55 
^ -0.60 0.80 J 



U2 



/ ^ 0.5 0.5 \ 

n 1 1 



(26) 



The first mixing matrix of Eq.26 corresponds to the solution of the solar neutrino anomaly based on the MSW 
mechanism with larger mixing angles (LMA) and the second corresponds to the vacuum oscillation solution 
for very low values of the squared mass difFerence(Just So) [37]. 

In order to investigate how the HEIDELBERG-MOSCOW experimental data constrain the parameter 
space of the supersymmetric seesaw model, we compute the lifetime of neutrinoless double beta decay in 
nuclei ^^Gg- In Fig-4, we plot the neutrinoless double beta decay lifetime of nuclei '^^Ge versus the lightest 
T-sneutrino mass m.g. For the CKM matrix of the lepton sector, we take Um = Ui. Fig. 4(a) corresponds 



to /i^ 



m 



- 1., = 10 GcV, and Fig. 4(b) to /13 = 0.1, = 10 GcV. The other relevant parameters are: 

4 X lO^GcV, tan/3 = 20 (solid lines) or tan/3 = 30 (dashcd-lincs). From Fig. 4, wc find that the 

"2 

theoretical prediction on the {PP)ov decay half life time of nuclei ^^Gg can meet the experimental data as long 
as the lightest r sneutrino has a mass less than ITeV. As the lightest r sneutrino mass increases to a certain 
value, the theoretically calculated values of the {(3(3)oi, decay half life time of nuclei ^^Gg are larger than the 
experimental upper limit. When we set the lepton CKM matrix as \J m = U2, and the other parameters to 
be the same as in Fig. 4, we re-calculate the relation of the half-life time versus ni- and the results are shown 
in Fig.5. The trend observed in Fig.5 is similar to that in Fig.4 and the similarity is understood as a property 
of the model. 



6 Discussion and Conclusion 

In this work, we analyze the neutrinoless double beta decay of nuclei ^^Gg in the supersymmetric seesaw 
model. With some assumptions on the model parameter space, we can naturally explain the solar and atmo- 
spheric neutrino experiments with the oscillation scenario and simultaneously the HEIDELBERG-MOSCOW 
experimental result of the ^^Gg neutrinoless double decay. In the numerical analysis, we assume that the 
masses of the three generation right handed neutrinos are degenerate, i.e. = m" = = m^, the as- 

sumption can simplify our calculation to a certain extent. Meanwhile, we take /i^ = 0, ^2 = To- choice 
corresponds to the hierarchical class solutions of neutrino masses in the three-generation models. As for the 
lepton CKM matrix, we choose two typical mixing matrices. In fact, we may have some other choices on 
the supersymmetric seesaw model parameter space to understand the neutrino oscillation together with the 
HEIDELBERG-MOSCOW neutrinoless double beta decay experiment. An important point is noted that for 
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each special assumption, we must have at least one sneutrino of middle energy scale, otherwise the data of 
the three experiments cannot be simultaneously fitted in this supersymmetric sea saw model. 

Being more explicit, in the preferred three-generation neutrino model, the contribution from the virtual 
Majorana-type neutrinos is suppressed by a factor < nii, > /pp of about 10~^^ and as well as the typical loop 
suppression factor at the electro-weak scale. Thus the SM with right-handed neutrinos cannot meet the data 
because of the unavoidable suppression factor. In contrast, as discussed in the text, the lightest sneutrino has 
a Majorana-type component which violates the lepton number and contributes to the neutrinoless double beta 
decay and it does not suffer from the strong suppression factor < rrii, > /pp- Generally, the fraction of this 
component is about 10~^. In addition to the typical loop suppression, the total suppression factor is about 
10~^ at the amplitude level. This value (or just the order of magnitude) can meet the data of the solar and 
atmospheric neutrino experiments and the HEIDELBERG-MOSCOW {(3(3)oi/ data simultaneously. Therefore 
the result implies that the observed data in the solar and atmospheric neutrino experiments can be explained 
by the oscillations mechanisms of the three-generation neutrinos, whereas the neutrinoless double beta decay 
is due to the contribution from the SUSY particles and the applied model is the SUSY see-saw model with 
the right-handed neutrinos. On other side, thus one can expect that the (PP)ou decay data greatly constrain 
the parameter space of the model, even though there still is large free room in the space which should be 
restricted by further experiments. 

Our conclusion is that the chosen parameter space of the supersymmetric model with right-handed neu- 
trinos can naturally explain the observed solar, atmospheric neutrino experiments and the HEIDELBERG- 
MOSCOW {PP)ou data. Even though the parameter space of the model still cannot be completely determined, 
the progress is noticeable. 

Introducing a "sterile" neutrino or a CPT-violation term in the neutrino sector, we can also accommodate 
all the lepton flavor violation processes [44]. How to distinguish between the supersymmetric seesaw model and 
those models is an interesting subject. Because the "sterile" neutrino participates in the weak interaction only 
via the mixing with the normal neutrinos, accurate measurements and systematic analysis on the r, /x rare 
decays may indicate the difference between the supersymmetry seesaw and "sterile" neutrino models. One 
plausible measurement which can distinguish the supersymmetry seesaw model from the CPT-violation model 
is the detection of neutrino magnetic moments in accurate experimental measurements (if possible). Because 
in the SUSY seesaw model, the neutrinos are purely Majorana-type which cannot have non-zero magnetic 
moments[45], by contraries, the neutrinos which reside in the CPT violation model can possess a Dirac-type 
neutrino component, thus can accommodate a non-zero magnetic moment. Obviously, such measurements 
would be extremely difficult, however, on the other side, almost all experiments on neutrinos are very difficult, 
so we lay our hope on the future developments of physics and technology. 
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A The sneutrino mixing under MIFF assumption 

With the MIFF assumption, the sneutrino mass matrix reduces to (one-generation) 

Pi Pi \ 

Pi Hi+hf^l + hKi 



with 



Defining a symbol 



m 



m 



(27) 



1 - 

w 

h^i^-^ + Km„iV2^A^^ (28) 

the three eigenvalues of sneutrino mass matrix are obtained as 

^li = I + ^hf^l + ml + + - 2A^) , 

li = \ {^ti + + + + 2^, + 2A^) , 

rr?=-(m'^^+2W\l + 'm?-BA. (29) 
Assuming the relations <C m? ^ + 2h^^V2 + or <C -B", Eq.29 is simplified as (up to order 

„J I nI I 

2 , 

=t,i 2~ ' 



1 



r^\=l(jnl+2hfvl + ml+B^y 
The corresponding unitary matrix which dragonflies the mass matrix is also simplified as Eq.23. 



B The expression for rji 

%i(Mw) = ^\ . Sn2 i^ud) U^nUop{zf+'^'^^ + ^(6+-)i)(z(3+'')i + zij^+<')^)zf Z^^ zf zf* 

2 ^ xlluXa 



^s^z^ +c^z^) X. ^n(x.-x„)' 
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VaM = 7^{Kd) ?7^n?7op(Zf+H.- +z(6+m)i)(2{3+o),- +_2(6+o)i)^nfc^p^ 

a=Li,Uj,Lk,Ki 
w ^ w 



(Z-'JZ]^* - -^Z'Jzf*) + mlm^_ {z'J*Z}! + ^Z'J*Z^) {z^J*Z\^ 
+^Z'fz'_') + 4m^_ m^m^, [z^Z^ - ^Z^JZ^) {z'fz}_' + ^^'f^'-')) , 

E 



Xa In a^a 2 



n {Xb - Xa) 



w w 

m 



2 1 '' - 1 

mpm n rv 

%(Mw) = - 8^^?^^ (K*,)^^en^eo^y (^f + + 

•7-2i*\ /' •3-2i*\ ^ c ^rlm* , •72m*\ -r2m* -3-1/ -'- -y-Am* ■y-2l\ 

j;^ In 

^ _ n(^6-^a) ' 

mpm 
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w w 

+ 52^+ ^ ^ _ „ n(x.-x„)' 

a=Ui,K. ,i/feA, Ail, 

mpm „ ^ 

W WW 

ig-^w^A. +^w2^ Ji^w^^ +^w2^ J (,2+ -—2+ J ^ 2^ _ ^ n(^6-^a) ' 

mpm n An / \ 2 

w w 

mpm „ ^ 
w w 



(-zlrs^ + z^r^*nY(z]'*-^zl^*)(z]!*-^z^J*) y 

w 

a=g,Ui,K. ,i^fc,Kj >t^m a^b 



w 



^ + 52 + + 92 ^ I n(^6-^a)' 

mpm _ 2 

%i(Mw) = T^-r^^iKd) UenUeo2fZ^^Z}!-*Zlj%Z(^+->+Zl^+->){Zl^+"> +^(6+o)Z*^ 

V + 52 + ^ ^ + 52 + A ^ + V2 ^ + ^ no- - n (^6 - ^g) ' 
%.(/^w) = 2^ 6 7c t/ent/eo2^'*2^'^l™*^i/*(^f+")'* 
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W W w 

f-- - - n ( 



mpm p 



mpm „ 



111 S^a 



_n On - - - Uixb-Xa) ' 
mpm^Q 2 

??j2lMwJ - 2 277r(72 V«<^J ^^otVep^^ ^- ^- 

w w 



n Of, - - - n (a^b - a;a) ' 
mpm- 4a, n2 



^ /,, \ f 9 ^Cig ( ^j.^ \ •T'li <rlfe ^li* •T'lfe* -T-ln* •T'li* /• •7(3+o)m* i <r(6+o)m*\ 

rikAt^w)- 2 qT,Vud) ^eoUep^^^^^^ ^- ^- + ) 

w 

^^(3+p)m* I ^(6+p)m*-|/^ln K2;2n\/-^ll 2;2l\ X^lnXa 



52 52 ^ ft - ff^- - - n(a;6-a^a) 



^ \ _ ^-P^g 4(3;, ^^^^ gr4i gr4fc ^rli* ^Ifc* ^rln* ^rli* / ^(3+o)m* I -7(6+0 

%2(./^w) - — 2— ^eot^ep^^ ^- ^- l^, + 



')m* ^ 



m^ 9tt 



52 52 ^_ft.^__ _n(a;6-a;a) 



mpm „ 

a 



^^i(/^w) = -^^^^^{v:d) u.oUevZ'^z^izfzfz':^*z^^{zf^'^-)^* + zf'^<^')^*) 

W WW 
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mpm n 
m2 37rC2 



)m* ^ 



E 



n {Xh - Xa) ' 

(V* U Z^^ Z^^* Z^^* Z^^* Z^^ ( Z^'^^'^^^* -I- z^'^'^°^^* ) ( Z^"^^^^^ _|_ ^C^+p)??!"! 



a=Di,K'-j,K.^, ,Ul,Um.,Kn 



In Xa 



i v*.Yu,„Ue„z^^z'^'z^"z^^*z^'''*zi^(z^^+"'>''''* + + z^^+p> 



mpm. 4as 



m4 Qtt 



E 



X^ InXa 



(31) 



Where Xj = The coefficient Ci{^^) {i = I, 2, • • • , 14) are defined as 



C2(/"w 
C3(/Ww 
C4(Mw 



^ai (/^w ) + ^bl (/^W ) + %2 (Mw ) + ^Cl (/^w ) + %2 (Mw ) + %2 (/^W ) + (/^W ) + %1 (/^W ) ' 

%2 ifJ-W ) + ^712 (Mw ) + (/^W ) + %1 (/^W ) ' 

2%i(/Xw) - %2(/^w) + ^(/^w) + %()"w) + Vm2{lJ'w) ' 

^ |-%2(/"w) ' 

W/ ' 

'7j2()"w) +%2(M 

'VfeiC/^w) + %2(/^w) - 'Vcil/^w) + 7^e2(/^w) ' 



Cio(/Xw) = 7%i(/^w) - %2(/iw) ' 
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Cii(Mw) 




Cl4(/^w) 



<^13Uw) 




(32) 
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Figure 2: The short-range contributions to ipp) in MSSMRN (Part one 
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Figure 3: The short-range contributions to 1/3/3 in MSSMRN (Part two) 
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Figure 4: The lifetime of neutrinoless double beta decay of nuclei Ge versus the lightest r-sneutrino mass 
with leptonic CKM matrix U^^ = Ui and {a)h^ = 1, = lO^^GeV, 771,3 = 4 x lO'^GeV; {h)^ = 0.1, = 

lO^^GeV, = 4 x 10^ GeV. The solid line corresponds to tan/3 = 20 and the dash line to tan/3 = 30; the 
other parameters are taken as in the text. 
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Figure 5: The lifetime of neutrinoless double beta decay of nuclei ^^Gg versus the lightest r-sneutrino mass 
with leptonic CKM matrix U^^ = U2 and {a)h^ = 1, = lO^^GeV, 771,3 = 4 x lO'^GeV; {h)^ = 0.1, = 

lO^^GeV, = 4 x 10^ GeV. The solid line corresponds to tan/3 = 20 and the dash line to tan/3 = 30; the 
other parameters are taken as in the text. 
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